Investigations of Ra + properties to test possibilities of new optical frequency standards 
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The present work tests the suitability of the narrow transitions 7s 2 Si/2 — > 6d 2 D 3 /2 an d 
7s 2 Si/2 — ► 6d 2 D 5 /2 in Ra + for optical frequency standard studies. Our calculations of the lifetimes 
of the metastable 6d states using the relativistic coupled-cluster theory suggest that they are suffi- 
ciently long for Ra + to be considered as a potential candidate for an atomic clock. This is further 
corroborated by our studies of the hyperfine interactions, dipole and quadrupole polarizabilities and 
quadrupole moments of the appropriate states of this system. 

PACS numbers: 



Accurate time and frequency measurement is crucial 
for the advance of science and technology in many fields. 
This leads to a number of searches to find candidates for 
optical frequency standards. The current frequency stan- 
dard is based on the ground state hyperfine transition in 
atomic cesium and has a quality factor (Q) of 10 15 
Atomic spectral lines with high Q are generally interest- 
ing for standards, however good control over systematic 
line shifts will be essential. As a result of the remarkable 
advances in the field of ion trapping and laser cooling, 
single ions like Hg+ , In+ [| , Ca+ fc| , Sr+ @, Yb+ Q, 
Cd+ 0] and Ba+ @] are particularly interesting as they 
can be localized using their electric charge rather than 
light forces, which is necessary for atom trapping. Very 
accurate measurements have been performed on Hg + and 
Sr + , where Q exceeds 10 17 . Some of the major systematic 
errors associated with the clock frequency are the Stark 
effect, Zeeman effect and quadrupole shifts due to stray 
electric fields in the ion trap [9( . These errors can be es- 
timated from high precision theoretical studies of hyper- 
fine structure constants, polarizabilities and quadrupole 
moments of the appropriate atomic states. Indeed, stud- 
ies of these quantities are also essential for parity non- 
conservation (PNC) studies 13, Hi- Some of the above 
mentioned errors can be eliminated by considering the 
clock transition between suitable hyperfine states |2j. 

An experiment is in progress at KVI to search for an 
suitable optical frequency standard by measuring the fre- 
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quency of either 7s S1/2 

6c? 2 D 5 /2 transitions in Ra + . A similar experiment is 
also being planned at IACS [12]. In this paper, we re- 
port our theoretical studies on the feasibility of these 
transitions for the optical frequency studies in Ra + . In 
the case of Ba + it has been pointed out that PNC and 
optical frequency standards experiments share many fea- 
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FIG. 1: (color online) Schematic diagram of energy levels of 
Ra + with transitions for possible optical frequency standards. 



tures in common The techniques used in the Ba + 
experiments can be extended to Ra + as the electronic 
structures of the two ions are similar. However, Ra + 
has one important advantage: the low lying transition 
wavelengths (see Fig. HJa)) of this ion are in the opti- 
cal regime making them more easily accessible than their 
counterparts in Ba + . Although, it appears that these 
wavelengths can be measured very precisely using mod- 
ern spectroscopic techniques, it is however necessary to 
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find out which transition is the most suitable for optical 
frequency standard. This can be decided by the experi- 
mentalists from the knowledge of different physical quan- 
tities of the involved states that can be used to control 
the sources of error. 

First of all, one must determine which of the isotopes 
of Ra + merit consideration for optical clock studies. In 
this context, it is worthwhile to note that only 223 Ra and 
225 Ra have half-lives of a few days (~ lOdays) and these 
isotopes are therefore obvious choices. However, they 
have different nuclear spins (Is); the former has 7 = 3/2 
whereas the latter has 1=1/2 and this results in different 
hyperfme splittings. One has to take into account the 
various systematic errors while considering both these 
isotopes. It is possible to eliminate the quadrupole Stark 
shift by considering the transition between the hypcrfine 
(F) states such as \6s(J = 1/2), I = 3/2; F = 2) -> 
\5d{J = 3/2), I = 3/2; F = 0) transition (see Fig. U£b)) 
for the frequency standard although knowledge of the 
hyperfme structure constants and the polarizabilities are 
still required for these experiments. It is necessary to 
study the hyperfme structure constants, lifetimes and a 
few other spectroscopic quantities for the 7s and 6c? states 
of this system in order to assess the suitability of the 
proposed clock transitions. 

Electron correlation and relativistic effects must be 
treated accurately for Ra + . Relativistic coupled-cluster 
(RCC) theory; a size-consistent, size-extensive and an all 
order perturbation method is well suited for this purpose 
fl3j |. It has been successfully applied to determine accu- 
rately certain grou nd and excited states properties of Sr + 
[3] and Ba+ [H]. We employ the same method in the 
present study to obtain accurate results for Ra + . The 
presence of the non-linear terms in this method makes it 
challenging to obtain the ground and excited state wave 
functions for a large system like Ra + . We had observed 
earlier that these effects are important for accurate stud- 
ies 16[ in other heavy systems. In order to obtain the 



wave functions for Ra + , we solve the RCC equations 
considering single, double and leading triple excitations 
(CCSD(T) method). This involves the determination of 
10 7 cluster amplitudes self-consistently. This is one of 
the largest computations to date for obtaining the wave 
functions of an atomic system. 

The starting point of our work is the relativistic gen- 
eralization of the valence universal coupled-cluster (CC) 
theory introduced by Mukherjee et al. [17| which was 
put later in a more compact form by Lindgren [H| . In 
this approach, the atomic wave function for a single 
valence (v) open-shell system is expressed as 



|*„> = e T {l + ^}|$„), 



(0.1) 



where is the reference state constructed out of the 
Dirac-Fock (DF) orbitals of the closed-shell system (|$o)) 
by appending the valence electron orbital. Here T and S v 



are the excitation operators from the core and valence- 
core sectors (for example see 14 , 1|| for the second quan- 
tization representations of these operators and equations 
to obtain their amplitudes). The single particle orbitals 
in the present calculations are linear combinations of 
Gaussian type functions flit ]. 

The transition matrix element of a hermitian operator 
(O) corresponding to the initial state |^) and the final 
state 1^/) can be expressed using the RCC method as 



(O) 



if 



($i|{l + gl}0{l + g/}|$/) 



(0.2) 



where we define O = e r1 Oe T and N v = ($„|e Tt e T 



Sle Ti e T S v \§ v ) for the valence electron v. We calculate 
the above expression using the procedure followed in the 
earlier works [3, 15, [l6j]. The expectation values are 
determined by considering the special condition i = f. 

Lifetimes of the 6d states: It is necessary to know the 
lifetimes of the 6d metastable states to understand how 
reliably the proposed experiments can be performed in 
that time period. The lifetimes (in second (s)) of these 
states can be determined from the inverse of the total 
transition probabilities (A) . The net transition probabil- 
ities (in s _1 ) of the 6d states are given by 
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1.11995 x 10 18 
Wf + 1)A 5 

2.69735 x 10 13 
Wf + 1)A 3 



cE2 



cMl 



(0.3) 

(0.4) 
(0.5) 



where Sf^i = |0/i| 2 and A (in A) are the transition 
line strength for the operator O (in atomic unit (au)) 
and wavelength, respectively. These quantities depend 
on both the transition amplitudes and wavelengths, and 
they can be calculated using a single ab initio method. 
However, we use experimental wavelengths [20j j to reduce 
the errors in the determination of the lifetimes. 

Since there are no experimental or theoretical predic- 
tions of the lifetimes of the 6d states, we calculate the E2 
transition amplitudes using both the length and veloc- 
ity gauges in order to assess the numerical accuracies of 
the results. These results are given in Table U along with 
the Ml transition amplitudes and the lifetimes of the 6d 
states. We have used the E2 amplitudes in the length 
gauge as it converges faster than the corresponding val- 
ues in the velocity gauge. The errors are estimated from 
the discrepancies of the results obtained with different 
choices of bases. 
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Transition Ofti Of-i, Lifetime 

states (au) (au) (s) 

/ — > i Length Velocity 



|64/ 2 ) -> |7s 1/2 > 14.87(7) 14.77(22) 0.0024(2) 0.893(4) 

|6d 5 / 2 > -<• |7si /2 ) 19.04(5) 19.87(1.0) 0.301(3) 
+ M5/2) -» |6d 3/2 ) 8.80(4) 10.5(2.5) 1.546(1) 0.297(4) 



State 


Odf 




e„ 


e 


64/2 


3.48 


-0.01 


-0.51 


2.90(2) 


64/2 


5.19 


-0.02 


-0.65 


4.45(9) 



TABLE II: Quadrupole moments of atomic states in au. 



TABLE I: Transition amplitudes (in au) due to Ml and E2 
transitions in both length and velocity gauges. Length gauge 
results of the E2 amplitudes along with Ml amplitudes are 
considered for the determination of lifetimes. 



Using the RCC method, we find that due to the en- 
hanced role of electron correlation, core polarization ef- 
fects in particular, the Ml transition amplitude for the 
|6d 3 / 2 ) — » |7si/2) transition is 0.0024(2)eao where the 
DF value is ~ 10~ 5 eao. From the calculated E2 ampli- 
tudes of this transition, we obtain the lifetime of the 6d 3 / 2 
state as 0.627(4)s. Inclusion of the above Ml transition 
probability changes its value to 0.893(4)s; this change is 
around 30% of the total result and this finding is different 
from our earlier studies on similar states of other alkaline 
earth metal ions j2^|. However, like the other d 5 / 2 states 
in those systems, the lifetime of the 6d 5 / 2 state reduces 
from 0.301s to 0.297s after including the contribution of 
the Ml transition probability in the 1 6rf 5 /2 ) ~* 16^3/2) 
transition. 

Quadrupole moments of the 6d states: In order to es- 
timate the error in the frequency of the clock transi- 
tion arising from quadratic Stark shifts, it is necessary 
to know the quadrupole moments of the relevant states. 
The quadrupole moment of a valence state (u) is given 
by 



e(v) = (# v \o^\9 v 



(0.6) 



where (9 E2 is the E2 transition operator. We divide the 
the above expression into three parts as follows 



G(v) = Qdf(v) + Qcv(v) + B v (v). 



(0.7) 



Here ®df, @ct and Q v are the DF, core- valence and va- 
lence electron correlation effects. In Table [Til we present 
these contributions for the 6d 3 / 2 and 6c? 5 / 2 states. In 
this table, the difference between the total RCC result 
and the sum of all the above three contributions is due to 
the normalization of the wave functions. The quadrupole 
moment of the 7s state is clearly zero as the quadrupole 
moment operator is of rank two. Therefore, we determine 
these quantities only for the 6d states. 

As given in Table HH the dominant contribution 
comes from ®df followed by 0„, which contains core- 
polarization and pair-correlation effects to all orders, 
make significant contributions as in Sr + 1J| and Ba + 



|15j | . We have followed the same procedure as in the life- 
time calculations to estimate errors in these results. 

Polarizabilities: We determine the dipole polarizabil- 
ities for 7s and 6d states and quadrupole polarizability 
of the 7s state for our study. The static (oq(J v )) and 
tensor dipole (ct\(J v )) polarizabilities for the valence v 
state with angular momentum J v are given by 



aj(t>) 

and 
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= 4, 



' 30^(2^-1X2^ + 1) 
(j„ + l)(2j„ + 3) 



E(- 1 ) J " +Jfc+1 



Jv 1 Jk 1 |(J„|£>|J fc )| 2 
1 J v 2 1 E v — Ek 



(0.9) 



respectively, where D is the El operator. Similarly, the 
static quadrupole polarizability {o^{v)) is given by 



al(v) = -4 



\{J v \Q^\J k )\ 
E v — Ek 



(0.10) 



We have used the sum-over-states approach and experi- 
mental energies to reduce the errors in the calculations; 
the calculated energies were used obtained from the RCC 
method where the experimental energies were not avail- 
able. 

We express generally the polarizabilities as 

a(v) = a DF {v) + a c (v) + a cv (v) + a v (v), (0.11) 

where each term is defined similar to the corresponding 
terms of the quadrupole moment expression given in Eq. 
(|0.7[) except for a c which is the pure core orbital contribu- 
tion. We calculate a v contributions from the calculated 
intermediate states using the RCC method. However, 
ct c and a cv are calculated using the second order many- 
body perturbation theory (MBPT(2)), where the residual 
Coulomb interaction and E1/E2 operators are treated as 
perturbation. All these results are tabulated in Table HTT1 

We have obtained up to 10s, lOp, lOd, 9/ and 9g low- 
lying states using the RCC method to calculate the above 
quantities. Contributions from other higher states are 
accounted for using MBPT(2). They are just given as 
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Statp 


CV£)_F 


"C 


ry 

"CD 


KX v 


CVj. 


CVS 




7si/ 2 

6*^3/2 
6rf 5 /2 


99.39 
996 39 

C/jjvJ .Oil 

1208.13 


14.57 
14.57 
14.57 


-1.54 
— 1.31 
-7.53 


-4.34 
39.05 
77.18 


-0.02 
0.03 
0.03 


107.86 
978.66 
1292.38 




6d 3 / 2 
6d 5 /2 


-220.67 
-168.53 


1.21 
1.21 


0.76 
1.27 


65.38 
16.78 


-0.01 
-0.01 


-153.12 
-149.28 


2 
Q 


7si/ 2 


2920.83 


56.57 


0.34 


-436.04 


5.50 


2547.20 



TABLE III: Dipole and quadrupole polarizabilities in au. 

tail contributions (at) in the table. Using the expression 

|<* 7s |0 E2 |v]/ fe }| 2 



a, 



0,6d 



(78) 



E 



E 6 



-,(0.12) 



we obtain a\(7s) — 1037(7)etQ along with the correspond- 
ing a c contribution. This is usually necessary for the 
lifetime measurements of the 6d states. 



7S1/2 



(id 



■3/2 



6f/: 



5/2 



Tpl/2 
8P1/2 
7P3/2 
8f 3 /2 
5/5/2 
6/5/2 
5/7/2 
6/7/2 



3.28 
0.04 
4.54 
0.50 



3.64 
0.07 
1.54 
0.15 
4.47 
0.86 



4.92 
0.40 
1.31 
0.21 
6.21 
1.08 



TABLE IV: Important reduced El matrix elements in au used 
to determine the dipole polarizabilities. 



In Table IIV| we present the important reduced El 
matrix elements which are used in the determination of 
dipole polarizabilities. These results are in reasonable 
agreement with those of Dzuba et al. which are calcu- 
lated using another many-body approach plj . 

Hyperfine structure constants: Studies of these con- 
stants are important to investigate the underlying physics 
of the wave functions in the nuclear region, especially to 
estimate the errors of the PNC matrix elements [26[ . The 
magnetic dipole (Ah) and electric quadrupole (Bh) hy- 
perfine structure constants of the valence v state with 
angular momentum J v are given by 



and 



A h (v) = ^(V v \TW\* v ) (0.13) 
B h (v)=2eQ N (^ v \T^\^ v ), (0.14) 



respectively. In the above expressions, hn, 9i and Qn are 
the nuclear magnetic moment, gyromagnetic ratio and 



quadrupole moment, respectively. Explicit expressions 
and the single particle matrix elements of TW and 
are given in 12311. We have used gj — 0.18067 1241 and 
Q N = 1.254 [II for 223 Ra and gi = -1.4676 [111 for 
225 Ra in these calculations. 



223 Ra+ 


7s 1/2 
A h 


6^3/2 

A h B h 


6^5/2 

A h B h 








RCC 


3567.26 


77.08 383.88 


-23.90 477.09 


Expt. 

225 Ra+ 


3404.0(1.9) 






RCC 


-28977.76 


-626.13 


194.15 


Expt. 


-27731(13) 







TABLE V: Hyperfine structure constants in MHz. 

The trends of the correlation effects in the hyper- 
fine interactions of the 7s and 6d states in the present 
system are similar to the corresponding states in Ba + 
[lE Si- We have found 23%, 31% and 181% correla- 
tion contributions with respect to the DF results of Ah 
in the 7s, 6c?3/2 and 6^5/2 states, respectively. The core- 
polarization (CP) effect in the 6d 5 / 2 state is very strong 
and its contribution is larger than the DF result. This 
gives rise to the unusual behavior of the electron corre- 
lation effects. 

Conclusion: We have successfully carried out accurate 
calculations of the lifetimes, polarizabilities, quadrupole 
moments and hyperfine structure constants in Ra + using 
the RCC theory. Our calculated values of the lifetimes of 
the 6d states which are 0.893s and 0.297s, respectively, 
suggest that Ra + could be a suitable candidate for an 
optical frequency standard with an accuracy better than 
10~ 18 . The results of the different properties that we 
have calculated can serve as benchmarks to guide exper- 
imentalists. On the other hand, precise measurements of 
these quantities can also be used to test our method of 
calculation. 
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